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Abstract. We describe K-mouflage models of modihed gravity using the effective field the¬ 
ory of dark energy. We show how the Lagrangian density K defining the K-mouflage models 
appears in the effective field theory framework, at both the exact fully nonlinear level and 
at the quadratic order of the effective action. We find that K-mouflage scenarios only gen¬ 
erate the operator at each order n. We also reverse engineer K-mouflage models by 

reconstructing the whole effective field theory, and the full cosmological behaviour, from two 
functions of the Jordan-frame scale factor in a tomographic manner. This parameterisation 
is directly related to the implementation of the K-mouflage screening mechanism: screening 
occurs when K' is large in a dense environment such as the deep matter and radiation eras. In 
this way, K-mouflage can be easily implemented as a calculable subclass of models described 
by the effective field theory of dark energy which could be probed by future surveys. 
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1 Introduction 


Future surveys such as Euclid [1] will actively look for deviations from the A-CDM paradigm. 
They could be resulting from dynamical properties of dark energy or induced effects of mod¬ 
ified gravity on the cosmological background evolution and/or the growth of structure at the 
linear level[2, 3]. Recently a new approach has been put forward which captures the main 
features of both approaches in an effective field theory description [4]. This effective field 
theory of dark energy allows one to describe models of dark energy and modified gravity 
using an operator expansion breaking Lorentz invariance explicitly[5-12]^. The coefficients of 
each operator are an unspecified function of time whose behaviour characterises each model 
at stake. Of course, this approach suffers from an enormous degeneracy, i.e. the appearance 
of many functions of time, with a multiplicity that grows with the perturbative order. This 
means that the applicability and the comparison with observational results of this approach 
is restricted to low orders of perturbation theory (typically once considers the background 
and linear levels). On the other hand, this provides a very general framework. For simple 
cases, such as the ones used in the Planck analysis [17], one can hope to extract information 
about dark energy and modified gravity scenarios from the comparison with data. 

The effective field theory approach also suffers from the fact that gravity acts on all 
scales, therefore implying that a perturbative approach at the cosmological level will usually 
not be able to unravel properties of dark energy or modified gravity models from the largest 
scales in the Universe to the Solar System ones. Such an extension is nevertheless crucial 
as one cannot be content with a large-scale description of dark energy or modified gravity 
which cannot be shown to be compatible with local measurements in the Solar system [18, 19]. 
The traditional approach to this issue, which involves building explicit models and screening 
mechanisms [20], complements the effective field theory description. For theories involving 
only second-order equations of motion, the compatibility between dark energy or modified 
gravity with local tests relies on the known screening mechanisms: chameleon[21-23] and 
Damour-Polyakov [24-28], K-mouflage [29, 30] and Vainshtein [31]. 

In this paper, we focus on the K-mouflage mechanism [32-37] and examine its description 
in terms of the effective field theory of dark energy. K-mouflage models are 7F-essence models 
[38, 39] complemented with a universal coupling of the scalar field to matter. The coupling to 
matter changes the cosmological dynamics of the models compared to iP-essence [32]. This 
can be exemplified by the looser link between the equation of state of dark energy rcde and the 
speed of scalar perturbations Cg [37] with no imaginary Cg when w^e < — 1- Because we start 
from a specific scenario, the K-mouflage modified-gravity model, which is fully defined at the 
nonlinear level, we obtain the exact nonlinear effective field theory action and equations of 
motion, ft turns out that the effective field theory action up to quadratic order involves the 
first two derivatives of the K-mouflage function K{x), where x is the reduced kinetic energy 
of the K-mouflage scalar whose presence affects gravity. In particular, the role of iF'(x), the 
first derivative of K{x), is made explicit inasmuch as it appears explicitly in the modification 
of Newton’s constant [33] -|- ^7)0 no and needs to be large in screened situations. 

In the cosmological context, screening occurs in the early Universe when densities are large 
and K' 1. On the contrary, when dark energy and modified gravity effects matter at small 
redshift, the function K' tends to one and the K-mouflage field behaves like a massless scalar 

^ Other methods for parameterising dark energy/ modified gravity have also been introduced in the last 
few years for models with second order equations of motion [13-16]. The comparison between the advantages 
of these different approaches is beyond the scope of the present paper. 
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coupled to matter with a strength (3. The coupling function A plays a crucial role in the 
cosmological evolution of Newton’s constant and is constrained by the local time dependence 
of Newton’s constant in the solar system providing a link between the cosmological and the 
local properties of K-mouflage. 

The K-mouflage models, and therefore their effective action, can be reverse engineered to 
accommodate the screening requirement if' 3> 1 in the early Universe (to ensure that the dark 
energy is subdominant at early times) and unscreening in sparse ones such as large cosmolog¬ 
ical scales in the recent past of the Universe. We show how by specifying two functions U{a) 
and A{a) at the background level, we can reconstruct the whole K-mouflage theories in the 
large-scale regime that applies to cosmological studies, down to cluster scales. The function 
U{a) is similar to the time variable of the background cosmology in the K-mouflage models 
expressed as a function of the scale factor U{a) ~ t{a). This new parameterization is very 
convenient as it allows us to normalize the models today and to directly set the magnitude 
and the behaviour of the deviation from the A-CDM cosmology through the function A(a). In 
the case of K-mouflage scenarios, the cosmological background and the growth of large-scale 
structures only probe the semi-axis X > 0 whereas small subgalactic scales, such as the Solar 
System, correspond to large negative x- Therefore, these models are a simple example of 
theories where the cosmological behaviour is decoupled from the small-scale behaviour^. 

The dark energy effective action for models like K-mouflage does not carry new infor¬ 
mation which cannot be extracted from the original non-linear Lagrangian. In that respect, 
the purpose of our paper is twofold. First we show that the complete K-mouflage theory 
can be defined not only by its Lagrangian but also by two functions of the cosmological 
scale factor. Given these two functions, it is straightforward to reconstruct the original La¬ 
grangian, to derive the cosmological background evolution and to deduce the effective dark 
energy effective theory at any order. This reconstruction can be easily used at second order 
in the effective theory, which is enough to analyse all the linear observables of the models 
and eventually to compare them with data. Secondly, the reconstruction mapping does not 
require any knowledge of the initial conditions in the early Universe and only requires the 
use of the cosmological parameters like flmO now. This is a drastic advantage compared with 
integrating the dynamics from initial conditions as the parameters of the model such as Ai, 
the cut-off scale of the model, do not have to be adjusted by a trial and error procedure. 
This is particularly important when implementing the comparison with data in a systematic 
way using a modified version of GAME for instance. The phenomenological work using the 
method expounded here is currently under way [40]. 

The paper is arranged as follows. In section 2, we define the K-mouflage models and recall 
the relevant equations of motion, at the background and linear levels. In section 3 we describe 
the effective field theory of K-mouflage. In section 4, we consider the cosmological dynamics 
at the background and perturbative levels and we check that we recover the equations of 
motion obtained from the usual K-mouflage action. In section 5, we show how the whole 
cosmological regime of K-mouflage models can be reconstructed from the knowledge of two 
functions, which provides a convenient parameterization of such scenarios. In section 6, we 
present numerical results for models behaving like cubic K-mouflage. We then conclude. In 
the appendices we provide details about the scalar field description of the quadratic effective 
action and about the causality of the theory. 

^ One may be able to infer the negative-y behaviour of K{x) from its positive-y behaviour by requiring 
analyticity, but this is not very practical as small deviations on the positive a:xis can be associated with large 
deviations on the negative axis. 
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2 K-mouflage 

2.1 Definition of the model 

We are interested in K-mouflage models, which we define by the action 


5 = 


J 


M, 


PI 


R + M^K{x) 




( 2 . 1 ) 

where is the Einstein-frame metric and g^y the Jordan-frame metric, given by the relation 

between Einstein- and Jordan-frame line elements 


ds^ = A'^{ip)ds^, g^u = A^{ip)g^y, 


( 2 . 2 ) 


when we use the same coordinates in both frames, and d.(y?) is the coupling function of the 
model. Throughout this paper, we denote Einstein-frame quantities by a tilde to distinguish 
them from Jordan-frame quantities. The first and second parts in the action (2.1) give the 
gravitational and scalar field contributions, which are defined in the Einstein frame. The 
K-mouflage Lagrangian of the scalar field ip involves a non-canonical kinetic term, defined by 
a dimensionless function K{x) with 


2M^ ’ 


(2.3) 


where A4‘^ is of the order of the dark energy scale today. The third and fourth parts in the 
action (2.1) give the matter and radiation contributions, which are expressed in the Jordan 
frame, and -i/ii are various matter fields. 

At the background level, we have the FLRW metrics 



ds^ = g^ydx^dx^ = d^(r)[—J- dx^j 

(2.4) 

and 

ds"^ = g^ydx^dx'^ = o^(t)[— dr^ J- dx^j. 

(2.5) 

where r 

is the conformal time, and the scale factors are related by 



a = A{(p)d. 

( 2 . 6 ) 


Throughout this paper, we denote background quantities with a bar. 


2.2 Dynamics of K-moufiage cosmologies 

We briefly recall here the equations of motion that govern the dynamics of K-mouflage cos¬ 
mologies in the Jordan frame. We refer the reader to [32-34, 37] for the derivations from the 
action ( 2 . 1 ). 

2.2.1 Cosmological background 

In the Jordan frame, the Friedmann equation takes the form 

- - d In A 

3i\TpjF7^(l — €2)^ = d + Prad + ptf with Mpi = A ^Afpp 62 = — -, (2-7) 

aino 
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where we have introduced the scalar field energy density given by 

^2 


p^ = A M (-K + 2xK ) with y = ( — 


dip'' 


( 2 . 8 ) 


Here t is the Jordan-frame cosmic time, with dt = adr, and H = dlna/dt the Jordan-frame 
Hubble expansion rate. The matter and radiation energy densities obey the usual conservation 
laws, which give 

(2.9) 


-/ \ PO A - / \ Prado 

p{a) = ^ and Prad(a) = 


a’ 

and the Klein-Gordon equation satisfied by the scalar field writes as 


d_ 

dt 


A-V 


dt 


= —a^p 


din A 
dp 


Ppo 

'Mpi’ 


where we defined 


/3(a) = Mpi 


din A 
dp 


( 2 . 10 ) 


( 2 . 11 ) 


2.2.2 Linear perturbations 

At the perturbative level, at linear order in the metric and density fluctuations, the Einstein 
equations can be analysed in the Newtonian gauge, 

ds"^ = a^ [-(1 + 2d>)dT2 + (1 - 2T)dx2] , (2.12) 

where the two Newtonian potentials $ and T are not equal. In the non-relativistic (<1> <C 1, 
'h <C 1 and ^2 <C 1) and small-scale {k/aH /§> 1) limits, the Einstein equations lead to 


$ = (1 -h ei)Tr 


2/32 V2 

T = (l-ei)TN, with ei = ^, —Tn = 

K' 


Sp 


(2.13) 


In the last equation in (2.13), we used the fact that fluctuations of the scalar field energy den¬ 
sity are negligible as compared with fluctuations of the matter density [34, 37], 5p^p <C 5p, in 
the small-scale limit. In a nutshell, this springs from the Klein-Gordon equation, which implies 
that 5p ~ /K'Mpik'^)6p when k/a S> H. This leads to dp^p ~ aP/K'k'^)5p <C dp. 

At the linear level, the effect of screening appears through the factor ei which measures 
the deviation of the two Newtonian potentials ‘h and T from Tn. When ei <C 1, gravity 
is not modified and the two Newtonian potentials coincide. It turns out that ei 7 ^ 0 is also 
responsible for the change in the growth of structure for K-mouflage models. Indeed the linear 


growing mode of the density contrast, D = ^, obeys the evolution equation 


d^D 


1 dH\ dD 


dt J dlna + 


+ 2 + ^ — 


(2.14) 


d(lna )2 

where the gravitational force is multiplied by the factor (Ij-ei) as compared with the A-CDM 
cosmology, i.e the effective Newtonian constant is 

^N,efr 




— 1 + Cl) 


(2.15) 


corresponding to a modification of gravity. At the linear level in perturbations, screening 
occurs when ei <C 1, which corresponds to backgrounds where K' ^ 1, as for the early- 
Universe cosmological background. It has to be noticed that the same factor ei appears in 
the static case and (2.15) is also valid. In this situation too, screening occurs where ei <C 1, 
which turns out to be inside the K-mouflage radius where GR is retrieved around a static 
object. 
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3 Effective field theory 


We now describe how the K-mouflage model defined by the action (2.1) can be recast in 
the effective field theory language. This can be useful to compare theoretical predictions 
with observations of the background and linear regimes, as data analysis may express such 
cosmological measurements in terms of the effective field theory functions in order to remain 
as general as possible. Then, within each class of scenario the parameters of the models must 
be determined to recover these few measured factors (if such a solution exists). 

One should be aware that the dark-energy effective-field-theory language which will be 
used throughout this paper is only a book-keeping device to analyse cosmological pertur¬ 
bations around a given background cosmology. In the case of K-mouflage, the background 
cosmology is defined by the non-linear K-mouflage evolution equations whilst the K-mouflage 
cosmological perturbations can be cast in the dark-energy effective-field-theory mould. As 
such the dark-energy effective field theory is only an expansion scheme for the cosmologi¬ 
cal perturbations of a classical theory, which is either explicitly known as in the case of the 
K-mouflage models or parameterised in a model-independent approach by unknown time- 
dependent coupling functions for each operator appearing in the effective Lagrangian. In the 
case of explicit models, the operators that appear are not necessarily all possible operators 
allowed by the generic cosmological symmetries (that are kept in a generic model-independent 
parameterisation) and we shall find that only one class of simple operators is generated by 
K-mouflage models. 

The dark-energy effective field theory should be contrasted with the quantum effective 
field theory obtained by including radiative corrections to the classical Lagrangian of a given 
model. The quantum effective field theory can be obtained by expanding in perturbation 
theory the classical Lagrangian around a given background configuration and then calculating 
loop corrections to the classical action. The quantum effective field theory of K-mouflage 
models will be analysed in a forthcoming publication. In a nutshell, the dimension-full scale 
which plays the role of an effective cut-off for K-mouflage models is K' In this case, it 

can be shown that the classical K-mouflage Lagrangian is not renormalised and that for scales 
larger than the inverse of the effective cut-off, a small scale compared to astrophysical ones, 
the finite correction terms to the classical action are all negligible. At these low energies, the 
particles of the standard model have all decoupled and it is legitimate to neglect the effect 
of the matter loops on the K-mouflage action. Hence it is perfectly legitimate to take as a 
starting point the classical and shift-symmetric K-mouflage action to develop the classical 
dark-energy effective-field theory. 

3.1 Unitary gauge 

In the general perturbed case, the scalar field (/^(x, r) depends on both time and position. 
However, we can make a change of coordinate to the uniform-(/j slicing or unitary gauge, 
defining for instance the coordinates and by 

X(„) = X, r(„) = <^"^[¥?(x,r)], (3.1) 

where the subscript (tt) refers to the unitary gauge. Thus we assume that ip can be used as a 
clock (which may not be valid on small nonlinear scales). Then, by definition the scalar field 
p = (p{Ti^u)) only depends on the time Thanks to the invariance under diffeomorphisms 
of General Relativity, the volume element transforms as d'^x^^—g = d^X{u)^—g{u)^ the Ricci 
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scalar as R = and the kinetic factor as 


2 


X 


= -X*{rin))9H = 


(3.2) 


We also move to the unitary gauge in the Jordan frame of the matter action (2.1) which 
becomes 




S = I df^xi 


(u) 9(u) 




(u) + X*9(u)) 9{u)fiu) SradiF'fii/i 9{u)fj.i/)i 

(3.3) 


with now 


9(u)fMu 9{u)pLv nnd A j4[y?(T^^^)]. 
Expanding in the perturbation of the metric, this gives 


(3.4) 


/' 


S= drx 


(u) \ 9{u) 


Mil 


( 00 )^ 


oo ^(n)._ JOO 


Kn)) 


n=0 


(3.5) 


+ 'S’m(V’i)^ 9{u)fiu) 4“ 1 A g(u)^v 


where we have defined = dAK/dx^{x) the derivative of order n of the kinetic function 
K at the background value x- This action can be written under the general form of the 
Einstein-frame effective-field-theory action of Gubitosi-Piazza-Vernizzi [4], as 


I 


S = d xJ-g^u) 


Ml 


PI 


^{u) 


W- 'l r(~ 'l^^OO I ^niRu)) , ^ QO . 

KRu)) - c{Ru))9(u) + 2 ^ {99(u)) 


n=2 


(3.6) 


+ g(u)iiv/f {T{u))) T *S'rad(T)jt/) 9[u)fiul f 


with 


/(t(„)) = A-^ A(r(,)) = -M^K + x*K'-fiu)i c(r(„)) = M\.K', 

n>2: M^i^^))=M\-x*rK^^l 

When the model is a simple canonical scalar field with a cosmological constant, 

^(x) = -l + X, (3.8) 

we find that 

A = 7W^ c = M^x*, n>2: = 0, (3.9) 

corresponding to both the cosmological constant and the kinetic energy of the scalar field. 
The inclusion of higher order terms in x in the K-mouflage Lagrangian will imply that both 
A and c will take more complex forms, with A becoming time dependent. 
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3.2 Jordan frame 

We can also write the action (3.3) in the Jordan frame, using the transformation laws 


9^J.u = f9iiu, R = f 


-1 


R-3ainf--g^^'^d^lnfd,lnf 


(3.10) 


In the case of the unitary gauge, this simplifies because / = /(t(„)) and the action (3.3) writes 
as 


S = 




'9(u) 




f dlnf\ 


J" 9(u)fiu) J“ Rra,d{,R^UT 9{u)^h 


Expanding in as in Eq.(3.5), we obtain an expression similar to Eq.(3.6), 


S = d x^-g(^u) 


A It i r(T ('^(m) ) 

R{u) ^{^{u))9(u) ^ 


n=2 


n\ 


(^5(n) 


00 ^n 


+ SYa(^'4’ii g[u)^v) J“ RrndiFfii/ j 9{u)fiu) ^ 


with 


(3.11) 


(3.12) 


M2i = /m2j, A = fA , c = /c-^M2i(^^^ 


, M^ = /2-m4 


(3.13) 


The expression (3.12) gives the Jordan-frame effective-field-theory action that corresponds to 
K-mouflage models, in the unitary gauge. In the usual effective field theory approach, the 
action is only defined as a perturbative expansion up to some low order, such as the quadratic 
level n = 2 if we consider the equations of motion up to linear order in the field fluctuations. 
This provides a general framework, as at each order one considers all possible operators, that 
applies to many different scenarios. On the other hand, the number of operators grows with 
the order and the theory is not specified beyond the truncation level. 

In this paper, since we focus on the K-mouflage models defined by the fully nonlinear 
explicit action (2.1), we have an explicit closed form (3.11) for the effective action and we 
obtain the expansion (3.12) at all orders. Moreover, it happens that the K-mouflage action 
(2.1) is simple enough to generate only a small subset of all possible operators at high orders. 
In fact, at each order we only have the single operator {6g^^)^. This is a key signature of 
K-mouflage scenarios amongst the class of theories coupled to matter in the context of dark 
energy effective field theories. If observations showed that some other operators cannot be 
set to zero, this would rule out all K-mouflage models. 


3.3 Equivalence with the scalar field description 

Although not transparent, the dynamics of the models defined by an effective action such 
as (3.12) are equivalent to a modification of General Relativity (GR) by a scalar degree of 
freedom. The scalar field which appears in the original definition of K-mouflage models can 
be unravelled by performing the changes [45] 

Rn) ^T + 7r(x, r), ^ df,{T 7r)du{T + Tr)g^'', 
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in the effective action that we have just obtained. We have the non-linear transformation 
ip{x,T) = (p{T -|- 7r(x, r)) and the perturbations of the scalar field are captured by vr in a 
non-linear way. At the linear level we have 

6if = ^TT. (3.15) 

dr 

The metric 51 ^,^ is the metric associated to the coordinates (x, r), where r is the background 
conformal time, and vr = 0 for the cosmological background. This “Stueckelberg trick” (3.14) 
restores the diffeomorphism invariance of the action, as an infinitesimal change of time can 
be absorbed through vrasr —tt ^ tt — 

In the appendix A we consider this change at the level of the quadratic effective action, 
as in the usual effective field theory approach. Here, because we actually know the action up 
to all orders since we start from the explicit nonlinear action (2.1) associated with K-mouflage 
models, we start from the full K-mouflage action (3.11) in the unitary gauge and change the 
time slicing. Both approaches are equivalent and lead to the same dynamical equations at 
the linear order of perturbation theory. 

After the transformation (3.14), the Jordan frame action (3.11), which was written in 
the unitary gauge where the time was specified by (3.1) explicitly, becomes 

5 = / TT)Mp^gf^^df,{T + Tr)d^{T + vr) 

+ f{T + Tt)M'^K [f~^{T + Tr)x*ir + TT)g^''d^{T + 'K)dy{T + vr)] j-hS'm -k S'rad- 

(3.16) 

This gives the full nonlinear K-mouflage action in terms of the scalar field tt in the general 
coordinate system (x, r). In the appendix A, we analyse the action at second order in vr, 
using the conformal Newtonian gauge for the metric g^y. 

Because the scalar field vr has been re-introduced through the unitary gauge effective 
action (3.11) and the Stueckelberg trick (3.14), which have a very general scope but are not 
specifically devised for K-mouflage models, the action (3.16) is significantly more complex 
than the original action (2.1). In particular, the scalar field vr is not identical to the original 
scalar field y?, though at linear level it corresponds to dip. 

The explicit expression (3.16) suggests that it may be difficult to extend the effective field 
theory analysis to high perturbative orders, unless one can perform some field redefinitions 
within specific classes of models, as even simple models such as K-mouflage give rise to many 
terms as one expand in vr. On the other hand, effective field theories are mainly used to 
study the linear regime or low orders of perturbation theory for which the number of terms 
is reasonable. 


rfln/ 

dr 


(r -k vr) 


4 Cosmological dynamics from the effective field theory actions 

We now describe how the equations of motions that determine the K-mouflage cosmology at 
the background and linear levels, which we recalled in section 2 . 2 , can be recovered from the 
effective field theory actions. 
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4.1 Background 


We can directly obtain the cosmological background dynamics from the unitary gauge action 
(3.12), using the metric (2.5) which also describes the unitary gauge in the homogeneous case. 
Because the Jordan-frame Planck mass Mpj depends on time, the functional derivative of the 
first term in Eq.(3.12) yields 


1 5{y/^fR) 




and the Einstein equations read as 


(4.1) 


M, 


PI 


- 2cJ°J° = T, 


(/ (^/+ 21^/) ] + (a - ^) 




-n, 

fiu 




(4.2) 


where and T^a-diiv are the matter and radiation energy-momentum tensors. Throughout 
this paper we denote with a dot derivatives with respect to the conformal time t, f = df /dr 
and / = d’^fjdr'^ [here we actually have / = df /dr]. 

The (00) component gives the Eriedmann equation 


3M|i 


(n^ + n 


dlnf \ 
dr ) 


0?p -\- O^Prad + O^A -|- C. 


(4.3) 


Changing time coordinate from the conformal time r to the Jordan-frame cosmic time t, 
defined as dt = adr (which is different from the Einstein-frame time t given by dt = ddr = 
A~^dt), this reads as 


3M|i (^H‘^ + =P + Prad + ^ + ^• (4-4) 

Then, substituting the expressions of A and c obtained in (3.13), we recover the Friedmann 
equation (2.7) recalled in section 2.2.1, which was derived in [37] in the Jordan frame from 
the original action (2.1). 


4.2 Perturbations 

To go beyond the background dynamics, we consider the fully nonlinear effective-field-theory 
action (3.16) expressed in terms of the scalar field vr. This allows us to derive the Einstein 
equations and the Klein-Gordon equation for vr. The Einstein equations read as 


Mv\[fGij.u — V^Vjy/ + g^y'0\f] — g^y 


-JMf.gP^dpT^d^T^ 7 + f M^K 


■\ 2 


(4.5) 


/ 


+ MpidpT^dyT^ I y 1 -2f M K a xd^r-^dnUT^ = + Ty^^dtiu, 


where we have defined the shifted time coordinate Ttt and the kinetic factor X, 

r^ = r + TT, X = -a^xg^’^^dpT^d^T^, 


(4.6) 









and we use f~^x* = In Eqs.(4.5) and (4.6), all explicit functions of time, /, /, and x, 
are evaluated at the shifted time which involves the scalar field vr. The kinetic functions 
K and K' are evaluated at the point X, which also involves the scalar field vr. 

The Klein Gordon equation for vr can be obtained similarly. 


Mpif 


R 3 


/ 


/ f 


— + ^ + T, 9 OyT-,, T “ ^9 


2 4 


/ 


X (X + 2nx'j g^''df,T.„duT^ = 2V/, 


/ f\ 

/ .\ 2 

/ 


+ 2//7W^K - fM^K'a^ 


I M - rM^K'a^xg>^''d,T, 


(4.7) 


where again the functions of time, /, /, a^, x nnd Xi nre evaluated at the shifted time Tjr, 
while the kinetic functions K and K' are evaluated at the point X. 

At the background level, the Einstein equations are obtained by putting tt = 0 and we 
retrieve Eq.(2.7). The Klein-Gordon equation reads when vr = 0 


SMli/ 


f2 f i 

- 1/2 I - 1 / J _ I 2 _ \_nj J 


+ 2ffa^M^K 


- \^ + mx + ^jx\- 2a^rM‘^K''xx = 0 , 




,2 f2 A /(4 


(4.8) 


which can be proved to be identically satisfied by using the Klein-Gordon and the Friedmann 
equations (2.10) and (2.7) of K-mouflage. Thus, the Klein-Gordon equation (2.10) of the 
K-mouflage field (/? in the Jordan frame and conformal time reduces to 


fM^K' 


3yX + 67^X + ^ 


+ 2fM^K"^^ = ^p, 


(4.9) 


where we substituted A = f The Friedman equation and its time derivative are given 

by 

= a^{p + prad + Pip) — 3Mpi/ y + 

and 


(4.10) 


- 2M|i/?7 = 


2 / P T Piad T Pp 


+ Prad +Pp] + ‘ 2 'M‘^^f ^ 


(4.11) 


where the scalar energy density and pressure are given by 

p^ = fM\2^K'-K), p^ = fM^k. 


(4.12) 


These identities combine to yield the Klein-Gordon equation (4.8) for vr at the background 
level. 
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At the perturbative level, to linear order in the equations of motion, the Einstein equa¬ 
tions can be analysed in the Newtonian gauge (2.12). The time-time component of Einstein’s 
equations (4.5) can be obtained by noting that 

Goo = 6Goo = -endr^ + 2VH, (4.13) 

where = VjV* is the flat-space Laplacian. In the sub-horizon limit of perturbation theory, 
the gradient terms dominate and we obtain in the quasi-static approximation 


V2 


a 


2 



2 M 2 j’ 


(4.14) 


which is nothing but the Poisson equation relating the Newtonian potential Tn to the matter 
perturbations (here we used that <C Sp in the small-scale limit). The spatial part of 
Einstein’s equations (4.5), using 


i/j: 5Gij = didji'if - 


(4.15) 


gives 


didi T - $ - 




= 0 , 


(4.16) 


in the non-relativistic limit <C 1, where v is the velocity of the matter fluid. At the linear 
order of perturbation theory and in the non-relativistic approximation we have therefore 


f f 

$ = 4'N-^vr, T = TN + ^7r, 


(4.17) 


where Tn is the Newtonian potential given by the last expression in Eq.(2.13). 

Similarly, the linearised Klein-Gordon equation in vr at first order can be deduced from 
Eq.(4.7) by using 

a^6R = (4.18) 


in the quasi-static and sub-horizon limit, leading to the linear relationship 

/ 4/32 

-vr = -^TN, 


(4.19) 


where /3 was defined in Eq.(2.11). This gives back Eq.(2.13), with a relative modification of 
the gravitational! potentials from Newtonian gravity given by the factor ei = 20^jK'. In 
particular, we retrieve that screening occurs when K' ^ 1. These results coincide with the 
ones in [37]. A comparison with the derivation from the second-order effective action is given 
in the appendix. This also leads to the relation between the effective field theory held vr and 
the original scalar held ip as 

^ -5^, (4.20) 


TT = — - 


/ Mpi 

where Sip = ip — (p \s the huctuation of the scalar held ip. Equation (4.20), which agrees with 
Eq.(3.I5), only holds at the linear level in the held huctuations. 
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5 Parameterizing K-mouflage 


The K-mouflage models are usually defined by the kinetic function K(x) and the coupling 
function A{(p) that enter the K-mouflage action (2.1). However, we may instead look for a 
parameterization of these K-mouflage models by two functions of the scale factor a, which 
would provide an implicit dehnition of K{x) and A{(p). This is similar to the tomographic 
approach developed in [41] for chameleon models where, instead of the function f{R) in 
/(i?)-theories or the scalar held potential V{^) and coupling A{lp) in Dilaton and Symmetron 
models, one dehnes the model through m^(a) and /3(a) [where m? = d‘^Ve?ildLp^ and /3/Mpi = 
d\i].A/dip\, which allows one to reconstruct V{(p) and A{(p). 

5.1 Properties of the K-mouflage fields and characteristic functions 

Before devising such an implicit dehnition, it is useful to recall a few properties satished by 
the kinetic and coupling functions. We refer the reader to [32-35] for the derivation of these 
constraints. Over all the cosmological evolution we require 

H > 0, /?' > 0, ^ < 0, /3 > 0, 62 < 0. (5.1) 

The condition H > 0 corresponds to the constraint that A does not vanish, so that the 
conformal transformation (2.2) is always well dehned. The constraint K' > 0 ensures that 
there are no ghosts. The scalar-held Lagrangian in the action (2.1) being even in y?, the sign of 
(f can be absorbed in the sign of (5 and we choose the convention /3 > 0. The coupling strength 
/3 does not vanish and keeps a positive sign if the scalar held is a monotonic function of time. 
In practice, the held excursion in units of the Planck mass is modest, from the primordial 
universe a —)• 0 until now a = 1, so that /3 does not vary much. Then, the combination K' > 0 
and /3 > 0 implies d(p/dt < 0 (as can be seen from the Klein-Gordon equation) and 62 < 0. 

In addition, we have the high-redshift behaviours 

a —^ 0 I A —y 62 —^ 0, K —y -l-oo, x —^ Too, ip —y ip^- (5.2) 

The hnite limit of A means that the Einstein and Jordan frames become equivalent at early 
time so that we can recover the standard cosmology, without a hfth force. This also implies 
that 62 goes to zero. The condition K' —)• Too ensures that the dark energy density becomes 
subdominant with respect to the matter density at early time. The marginal case where K' 
goes to a constant K* 3 > 1 , so that in the early matter era the dark energy density is a small 
hxed fraction of the matter density, is also possible but we do not explicitly consider this 
case here (it can be obtained as a limit of the range of models we consider). The condition 
X Too is a natural consequence of K' —)• Too if we consider kinetic functions that have no 
singularity at a finite value of x- Note that because we look for a tomographic reconstruction 
of the K-mouflage model, we restrict ourselves to models where x{a) is a monotonic decreasing 
function of a to avoid reconstructing ill-dehned multi-valued functions. Although the time 
derivative of the scalar held, d(p/dt, diverges at early time like the scalar held goes to a 
finite value. This is again due to the requirement that we recover the standard cosmology at 
early time, which implies that the scalar kinetic energy cannot grow too fast at high redshift. 

At late time in the future, far in the dark energy era, the scalar field ip will again converge 
to a hnite value and we have 

a—>- 00 : Too; A —>■ Aooi 62 —>■ 0. (5-3) 
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Finally, the kinetic function must satisfy the property 

\/^K'{x) is a monotonic increasing function of x over X > Oj K' + 2xK" > 0, 

r- (5-4) 

and y'x-^^(x) +00 for x ^+oo- 

This ensures that the system is well defined up to arbitrarily high redshifts, when p — )■ +oo, 
and that there are neither small-scale instabilities nor causality problems. 

In this paper we focus on the cosmological background and large-scale perturbations, 
which correspond to x > 0. The kinetic function must also satisfy similar properties on the 
semi-axis x < 0; which corresponds to the small-scale quasi-static regime, but we do not need 
to consider this here. 


5.2 Building a simple parameterization in terms of the scale factor 
5.2.1 The parameterisation 

In order to build an implicit parameterization we must consider the equations of motion 
satisfied by the system at the background level, namely the Friedmann equation (2.7), the 
usual conservation equations for the matter and radiation components (2.9), and the Klein- 
Gordon equation (2.10) satisfied by the scalar field. From the second relation (2.8) and the 
sign convention (5.1) we have 


dip 

dt 




dip 




A din a AH 

and the Klein-Gordon equation (2.10) reads 

d 


din a 


A-^a^^2xM^K' 


J_ ^ e2AH 

Mpi 


€ 2^/00 


(5.5) 


(5.6) 


This equation suggests that a simple parameterization of the K-mouflage model is provided 
by the two functions of the scale factor 

U{a) = a^\/^K' and A{a), with AQ = A{a = l) = \, (5.7) 

where we choose to normalize the conformal factor A to unity today. Given A{a) we obtain 
at once 62 ( 0 ) from its definition (2.7) and the Klein-Gordon equation (5.6) yields 




PO _ €2^^ 

2t/(-362 + S) 


(5.8) 


This gives an explicit parameterization of x(a). 


5.2.2 Normalising Ai 

The full reconstruction of the Lagrangian K{x) the coupling function A{ip) can only 
be achieved when the scale Ai has been specified in terms of the present-time cosmological 
parameters. This can be achieved using the Friedmann equation (2.7) which writes as 


772 

-”0 


^2 

(1 - 62)2 


17 


mO 


+ 


17 


radO 


+ 17(po^ 


K - 2^K' - 
Ko - 2|o^^J ’ 


(5.9) 
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where we defined the cosmological density parameters as Oj(a) = Pi(a)/pcrit(fl) and /3crit(o) = 
3Mp^{a)H‘^{a). In particular, from the scalar field energy density (2.8) we have 

- - - 

^mo (—-K'o + 2xo-I^o)- (5.10) 

Po 

To close the system we need to specify the scale in terms of the cosmological param¬ 
eters today, as we wish to normalize the K-mouflage model at z = 0. To make the Friedmann 
equation (2.7) look closer to the A-CDM cosmology, we introduce an effective dark energy 
density defined by 

3MpiA = P + Prad + Pde- (5-11) 

(However we keep the explicit time dependence of the Planck mass on the left hand side as 
the Planck scale is intrinsically time-dependent in the Jordan frame and this effect could be 
detected by Solar System experiments.) This yields explicitly 


Pde — Pip T 


2^2 - ej 
(1 - ^ 2)2 


(p T Prad T Pp) ■ 


The Jordan-frame dark-energy density evolves as 


= -3i7 (pde + Pde), 

where we have defined the Jordan-frame dark-energy pressure as 


(5.12) 


(5.13) 


Pde =Pp + 


(Prad T Pp) T 

1 — 62 



2 de2 \ P + Prad + Pp 
1 —e2dlnay 3(1 — 62 )^ 


(5.14) 


Notice that the coupling to matter modifies the usual expressions for the dark-energy pressure 
and energy density from the K-essence models. Similarly, the effective dark-energy equation 
of state in the Jordan frame is simply defined as 


'^de — Pde/ Pde- 


(5.15) 


This differs significantly from the AT-essence equation of state w^p defined as 


= ^ 

2lK'-K 


(5.16) 


In particular, the speed of scalar perturbations in K-mouflage models is not directly related 
to (1 J- tCde) like for K'-essence models [43]. We will see that models can cross the phantom 
divide Wde < — 1 in the past of the Universe without any instability, as remains positive. 
We can use the Friedmann equation to normalize 


Um T ^rad T ^de — 


Um T ^rad T 

(1 - 62)2 


(5.17) 


whence 


^(/jO — (1 “ 62 ,o)^^deO “ (262,0 “ f 2 ,o)(^mO + S^rado) with 62,0 = 62(0 — !)• (5.18) 
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This provides the value of the cosmological parameter today in terms of the standard cos¬ 
mological parameters {rimO) f^radO) f^deo} that would be derived from the Friedmann equation 
today assuming a A-CDM cosmology. In realistic models, which remain close to the A-CDM 
scenario and must satisfy Solar System tests, we have |e 2 | < 0.01 and fi^o — ^deo- This 
specihes the constant from Eq.(5.10), noticing that XqKq = \/^oUq from the definition 
(5.7) and using Eq.(5.8). Choosing the normalization Kq = —1 [as seen from the scalar-held 
Lagrangian in the action (2.1), we must choose a normalization of the kinetic function K to 
specify the scale this yields 


Po 


rimo 


^ 2,0 


Q, , dlnU 


a=l 


with Kq = —1, 


(5.19) 


and we recover the fact that Ad^ ~ pdeO- 


5.2.3 The reconstruction mapping 

Once the energy scale Ad^ has been obtained from the cosmological parameters, using Eqs.(5.18) 
and (5.19), and the values of the given functions e 2 and dlnW/dlna today, we directly obtain 
x(a) from Eq.(5.8). Next, from x(a) and the dehnition of 17(a) in Eq.(5.7) we obtain K'(a), 


K'(a) 


U{a) 

a^Vxia) 


This also gives K(a) from the integration of 


(5.20) 


K{a) 



U(a) dx 

a? v^x(o) 


(a) din a. 


(5.21) 


with the initial condition K(a = 1) = Kq = — 1. Inverting the mapping a —)• x(a), we can 
reconstruct K(x) from U{a) and A(a). 

Finally, the Friedmann equation (5.9) provides H(a), as a function of the Hubble rate 
today Hq. The value of the scalar held is obtained from the integration of the second 
expression in (5.5), with the initial condition (^o = 0 (which has no physical impact) 


= “X -to-***’'”' 


(5.22) 


The curves {x{a), K(a)} and {(p(a), A(a)} provide a reconstruction of the kinetic and coupling 
functions K(x) and A(^p), over the range probed by the cosmological background. This is 
a tomographic reconstruction mapping as the whole theory can be reconstructed from the 
knowledge of U{a) and A(a) as functions of the scale factor. 


5.3 Properties of the characteristic functions U(a) and A(a) 

5.3.1 Cosmological evolution 

The functions U(a) and A(a) that dehne the K-mouhage model through the parameterization 
described in the previous section cannot be fully arbitrary. They must obey several constraints 
that are associated with the properties (5.1)-(5.4). 

First, the Klein-Gordon equation (2.10) can actually be integrated as [32, 37] 

A-^U(a) = A-^a^^K'= [ dt^^^=, (5.23) 

Jo MpyV^ 
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that is, the integration constant must vanish in order to obtain a realistic early-time cosmol¬ 
ogy. Indeed, if the integration constant C of the Klein-Gordon equation does not vanish we 
obtain x ~ f-4/(2m-i)^ where K{x) ~ x”* with m > 1 in the highly non-linear regime, and 
~ ^ a Scalar energy density that grows faster in the past than the 

matter density. Thus, recovering the A-CDM model at early times implies that C = 0 [32]. 
Then, at early time we have d(p/dt —)■ —oo but (p goes to a finite limit, which implies that 
A and /3 also go to constants (this also means that the Einstein and Jordan frames become 
equivalent at the background level at high redshifts), and we must have U{a) oc t. This gives 

a <C 1 : U{a) (X in the radiation era and U{a) oc in the early matter era. (5.24) 

Next, since we focus on simple models where the various fields are monotonic functions of 
time, hence X is a monotonic decreasing function of a, the condition (5.4) implies 

U{a)/a^ is a monotonic decreasing function of o. (5.25) 

The condition U{a)/a^ —>• J-oo for a —)• 0 is already implied by the constraint (5.24). From 
the dehnition (5.7) and the relation (5.6), using €2 < 0, we obtain 

U > 0, U/A? \s a. strictly increasing function of a. (5.26) 

The function A(a), and its logarithmic derivative 62 ( 0 ), must obey the constraints (5.1) 
and (5.2). At early times we obtain from Eq.(5.8) ~ —e 2 /U. Since we look for a 

monotonic decreasing function x{'^)i with x J-oo for a —)• 0 , we obtain 

a <C 1 : —^ 2 /U is a decreasing function of a, —e 2 /U —t +00 for a —)• 0. (5.27) 

This means that €2 cannot go to zero too fast at high redshift. From (5.24) we can see that 62 
converges to zero more slowly than in the early matter era and than a? in the radiation 
era. 

At late times, the K-mouflage scenario leads to a A-CDM like accelerated expansion, 
driven by a constant dark energy density. The scalar held (p and the conformal factor A 
converge to constant values, t ~ ln(a) and U{a) ~ ln(a). Therefore, far in the dark energy 
era, we have 

a^l : f7(a) oc ln(a). (5.28) 

This asymptotic behaviour applies to the future; the current epoch, where dark energy is only 
starting to dominate the Universe, only sees the beginning of this logarithmic dependence. 
However, this means that over the range 0 < a < 1 the function U(a) should hrst grow 
as a^, and then as from (5.24), and further slow down at a < 1 to eventually reach 

the logarithmic growth of (5.28). In particular, choosing functions U{a) that are strictly 
increasing functions of the scale factor automatically ensures that the condition (5.26) is 
satisfied, because A{a) is monotonically decreasing (€2 < 0 ). For conformal factors A and 
strength /3(o) that do not vary much over the cosmological evolution, Eq.(5.23) actually shows 
that U{a) grows as U (a) ~ t(a) where t(a) is the cosmic time, and this should be close to the 
A-CDM expansion. 
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5.3.2 Normalising U(a) 

As shown in section 5.2, the parameterization of the K-mouflage model in terms of U{a) 
and A{a) allows us to normalize the cosmology to the present cosmological parameters, 
{rijnO) ^radO) ^deO)In addition, the function A is normalized to unity today and goes 
to a finite value A* > 1 at high redshift, with |e 2 | < 0.01 to satisfy Solar System constraints, 

A(0) = A*>1, A(l) = l, e 2 < 0 , |e 2 |< 0 . 01 . (5.29) 

This also means that A{a) is always close to unity, with typically (A* — A) < 0.1. 

When we define the K-mouflage model from its kinetic function K{x), as in the action 
( 2 . 1 ), we usually impose the low-y expansion K{x) = —1 + X + where the dots stand for 
higher-order terms, which sets the normalization of the model without any loss of generality. 
The factor —1 defines the normalization of the energy scale parameter The normalization 
to unity of the linear term defines the normalization of the scalar field (/? and of the kinetic 
factor X- We have already performed the first normalization in Eq.(5.19) in the case of the 
{17(a), A(a)} parameterization, except that we choose K[a = 1) = —1 instead of K{a —)• 
-|-oo) = —1. Thus, we still have one degeneracy associated with the second normalization, 
which we have not performed yet. Indeed, going through the reconstruction method described 
in section 5.2, we can check that the cosmology and our relations are invariant through the 
transformation U —)• XU, X ^ x/i 77' —)• X^K', K ^ K, ip ^ ‘^/'X- Therefore, we need to 
choose the normalization of 17(a) today, to fully characterize the model. Using this symmetry 
we normalize 17(a) by 


Kq = 1 , whence Uq = 


-Po^2,o 


M^2(-3e2,o + 


dlnU 

din a 


a=l)’ 


(5.30) 


where the second equation comes from the combination of Eqs.(5.7) and (5.8). Thus, the 
normalization of the kinetic function K (x) is not exactly identical to the one used in previous 
works [32, 34, 36, 37]. Instead of 77(x = 0 ) = —1 and K'{x = 0 ) = 1 we now have 
K{x = Xo) = cind K'[x = Xo) = 1- However, in practice we have xo Ij as can be seen 
from Eq.(5.23) which yields xo ~ and we must have /3 < 0.1 to satisfy Solar System and 
cosmological constraints [35]. Therefore, these two different normalizations are usually rather 
close. 

Finally, the standard A-CDM cosmology is recovered in the limit 62 —)• 0, as 17 
with the normalization choice (5.30), x ~ \^ 2 \, K' ~ 1, |i7 —i7o| ~ |e 2 [, (p ~ /3 ~ 

and IA — 11 ~ 1621 . 


5.4 What does 17(a) represent? 

It is interesting to see more precisely how the shape of the functions U (a) and A (a) are related 
to the shape of the kinetic function K{x). For a power-law large-x behaviour, 

X S> 1 : 77(x) ~ x”* with m > 1, (5.31) 

we have in the early matter era 

early matter era: ^ ^ a-^A‘2m-i) ^ U ^ ^ t, £3 ~ (5.32) 

and in the radiation era 

radiation era: ^ ^ a-^A'^rn-i) ^ U t, £3 ~( 5 , 33 ) 
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We can check that these behaviours satisfy the constraint (5.27) over the full range 1 < 
m < + 00 . We recover the behaviours (5.24) and we find that the function U{a) is largely 
independent of the kinetic function K{x)- In agreement with Eq.(5.23), it is mostly related 
to the expansion history of the universe as 

U{a) ~ t{a) (5.34) 

and should not deviate much from the “A-CDM shape” 17A-CDM(fl) that is obtained by setting 
A and /3 to constants in Eq.(5.23), unless one considers coupling functions A{(p) with steep 
slopes. The main freedom, associated with the shape of K{x), is contained in the function 
€ 2 ( 0 ). In particular, its power-law exponent in the scale factor, €2 ~ is directly related 

to the exponent of the kinetic function, with = 3(m — 1)/(2m — 1) in the early matter era. 


5.5 Factors ei and €2 

In addition to the factor 62 defined in (2.7), a key factor that enters the growth of large-scale 
structures is the function ei defined in (2.13), which enters the evolution equation (2.14) of 
the linear growing modes of the density contrast and describes the amplification of gravity in 
the unscreened regime. 

The function 62 ( 0 ) is directly obtained from the function A{a) within the parameteriza¬ 
tion {U (a), A(a)}, as €2 = din A/dln a. It is interesting to see how the characteristic function 
ei(a) is related to this parameterization. From the last Eq.(5.5) and the definition (2.13) of 
ei, we have ei = e^A"^xK'. On the other hand, from the definition (5.7) we have 
XK' = ■s/^Uja?' and from Eq.(5.8) we obtain x^' = —poe2A^/2A4^a^(—3e2 + d\n.U/d\n.a). 
Combining these two results gives 


€l(a) = - 62 


2 -362 + 


dlnC/l 
d\na ) 




(5.35) 


In the matter dominated era we have Ilm ~ 1 and dlnU/dlna ~ 3/2, see (5.24). This yields 


matter era: 61 ~ — 62 . 


(5.36) 


At late time in the dark energy era the ratio 61 /I 62 I grows because decreases in the 
denominator of (5.35). Therefore, in the K-mouflage models the functions 61 and 62 are equal 
for most of the cosmological evolution and measure both the drift of Newton’s constant (or 
the Jordan-frame Planck mass) and the amplification of the gravitational force in the growth 
of large-scale structures. 


6 Numerical results 

We now present our numerical results for several K-mouflage models defined by the parame¬ 
terization {U{a),A{a)} described in section 5. 

6.1 Explicit models 

In the following sections, we describe our reconstruction procedure and our analysis with 
a simple choice for the functions U{a) and A{a), with a few parameters. We shall consider 
three sets of parameters that cover different behaviours for the underlying K-mouflage models, 
within the constraints that must satisfied. The precise forms of the functions (6.1) and (6.4) 
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have no particular theoretical meaning. They are only the simplest choices that satisfy the 
previous requirements and are sufficient for our illustration purposes. 

For the function A{a) we take 


7 A > 0 : A{a) 


1 + aA — oiA 


{'jA + l)a 




7A + a 


( 6 . 1 ) 


This is the simplest function that goes to constants at both early and late times, a —?■ 0 and 
a —>■ oo, with a power-law behaviour at early times. The parameters 7^1 and va are 
related to the more physical quantities 62,0 and m by 


^ 2,0 — 


aAlAVA 
7A + 1 ’ 


VA 


3{m — 1) 

- with m > 1 , 

2m — 1 


( 6 . 2 ) 


where m is the exponent of the kinetic function in the matter era regime, defined by i^T ~ y*”, 
see (5.32). If we require the kinetic function to keep the same exponent in x ia th® radiation 
era we should change ua for a < Oeq, following (5.33). However, because the dark energy 
component is negligible in the radiation era and cannot be measured at these high redshifts 
(for the models that we consider here) we keep the same form ( 6 . 1 ) with a single parameter 
UA given by Eq.(6.2) at all redshifts. Thus, the parameterization (6.1) is determined by the 
three parameters {m,e 2 , 0 ) 7 A}) ciA being obtained from Eq.(6.2). As explained above, m is 
the logarithmic slope of the kinetic function at large x, 62,0 = din A/din a at a = 1 today, 
and 7 yi ~ 1 describes the transition to the dark energy epoch, where A converges to a finite 
value as in (5.3). In particular, we have 


A{a = 0) = 1 -|- a A, A{a = 1 ) = 1 , A{a = 00 ) = 1 -|- ua — aA( 7 A + l-Y^ ■ (6.3) 


The parameter € 2,0 measures the deviation of A from unity and the departure of the K- 
mouflage model from the A-CDM cosmology. 

For the function U{a) we consider the shape 


«(/ > 0 , 7 [/ > 1 : 


U{a) oc 


( ^/d) ln( 7;7 -I- a) au' 


(6.4) 


which comes from the simple combination 


lA-lu + a)' 

where Oeq is the scale factor at equality between the matter and radiation components. This 
is the simplest function that follows the evolution of t{a) through the radiation, matter and 
dark-energy eras, each term successively dominating in Eq.(6.5). In agreement with the 
requirements (5.24) and (5.28), it gives U{a) ~ first, followed by U{a) ~ and next 
U(a) ~ In a. The transition between the radiation and matter eras is set by Oeq whereas the 
two parameters ajj and xu describe the transition to the dark energy era and are of order 
unity. Finally, the normalization of U{a) is obtained from Eq.(5.30). 

In the following, we present our numerical results for the three K-mouflage models 
defined by the parameters given in Table 1. They have 62 = —0.01, which ensures that Solar 
system tests and cosmological bounds are satisfied, while other parameters are of order unity. 
We keep xu a-nd m unchanged, because the results are less sensitive to these parameters, 
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model 

au 

lu 

m 

£2,0 

7A 

1 

0.2 

1 

3 

-0.01 

0.2 

2 

0.2 

1 

3 

-0.01 

0.5 

3 

0.4 

1 

3 

-0.01 

2 


Table 1. Parameters for the three K-mouflage models studied in section 6. 



Figure 1. Left panel: background conformal factor A{a) as a function of the scale factor, for the 
three K-mouflage models defined in Table 1 where ja decreases from bottom to top. Right panel: 
background scalar field in units of the Einstein-frame Planck mass, (p/Mp\, as a function of the scale 
factor. 

and we vary ja and ajj. The main parameter that determines the cosmological behaviours 
is 'jA- Since m = 3, all three models satisfy a cubic form, K ~ in the highly non¬ 
linear regime. However, as we shall find below, these three sets of parameters correspond 
to significantly different kinetic functions K{x) and allow us to explore a wide range of K- 
mouflage behaviours. 

For the cosmological parameters we choose the Planck mission results, Hmo = 0.3175, 
f^deo = 0.6825 and Hq = 67.11 km/s/Mpc. 

6.2 Reconstruction of the kinetic and coupling functions 
6.2.1 Evolution with time 

From the functions U{a) and A{a) and the cosmological parameters today, we can compute 
the expansion history of the Universe, as described in section 5.2. We show our results for 
the three K-mouflage models defined in Table 1 in Figs. 1-3. 

We can see in the left panel of Fig. 1 that as 7 a decreases, from models 3 to 1, the 
deviation of the coupling function A from unity grows, in agreement with Eq.( 6 . 1 ) and A{a = 
0) = 1 -|- UA- This will lead to a greater deviation from the A-CDM cosmology, as the fifth 
force arises from the variations of the conformal coupling factor A. 

As seen in the right panel of Fig. 1, a greater cosmological range for A{a) also leads 
to a greater excursion for the field (p. Indeed, from the definition of /3 and 62 we have 
dip/dlna = Mpie 2 // 3 . Therefore, if /3 does not vary too much, a greater range of A, whence 
a greater |e 2 | at intermediate redshifts, leads to a larger time derivative of the scalar field. 
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Figure 2. Left panel: coupling strength /3(a) = M-p\dh\A/dif. Right panel: functions ei(a) = 
2fdf‘jK' (upper solid curves) and € 2 ( 0 ) = d\nA/d\na (lower dashed curves). 


Thus, as could be expected, a greater range for the coupling function A{(p) corresponds to a 
greater range for the scalar field (p itself. 

At low z, both A{a) and (p{a) converge for all three models. This is because we normalize 
all models to the same cosmological parameters and to the same slope € 2,0 = —0.01 of A{a) 
today. 

As seen in the left panel of Fig. 2, /3 remains about 0.1 for all three models and at all 
redshifts. Indeed, we typically have |e 2 I ~ ei = 2/3^/iF', see Eqs.(5.35) and (5.36). At late 
times we also have X — 0 and AT' ~ 1 , which gives |e 2 | ~ ei ~ 2/3^. Therefore, |e 2 ,o| = 0.01 also 
implies /3o ~ 0.1. Then, /3 remains almost constant with redshift because we chose a function 
U{a) that behaves like t{a), see Eqs.(5.23) and (5.24). This ensures a coupling function A[ip) 
that is smooth and “natural”, without sharp knees or oscillations that would introduce new 
scales. 

We show the two factors ei and 62 in the right panel of Fig. 2. In agreement with 
Eqs.(5.35) and (5.36), ei ~ —62 at high z and ei/|e 2 | grows at low z as l/Hm- The factors /3, 
ei and 62 converge at low z for all three models, because we normalize the models at z = 0 . 
For the model 1, which shows the greater deviation of A from unity at high z, the amplitude 
of 62 slightly grows from z = 0 to z ~ 4. At early time, a —)• 0, ei and 62 converge to zero for 
all models. For €2 this is ensured by our choice of the coupling function A{a), with €2 ~ —a'^^ 
from Eq.(6.1). Since at early times we have ei ~ —€2 this also implies ei ~ From the 
definition of ei in (2.13) this also means that K' —>• +00 at high z. The fact that ei is of 
order 1% implies that deviations of the linear power spectrum from the A-CDM prediction 
will be of the order of a few percents, see Eq.(2.14). 

Finally, we show the kinetic terms x and K' in Fig. 3. Again, both x and K' converge 
at low 2 ; for all three models. We can check that x and K' increase at high z. When K' ^ 1 
the K-mouflage mechanism comes into play with the dominance of non-linear effects in the K- 
mouflage Lagrangian. This cosmological K-mouflage screening ensures that the dark energy 
component becomes subdominant as compared with the matter density. However, in contrast 
with the constant dark energy density shown by the A-CDM cosmology, pde grows with z. 
We find that the smaller values of 7 ^ 4 , associated with larger A and (p, correspond to greater 
X and smaller K'. Indeed, at fixed U{a) we can see from Eq.(5.7) that K' so that 
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Figure 3. Left panel: background kinetic factor x- Right panel: first derivative of the kinetic 
function, K' = dK/dx- 




cp / 


X 


Figure 4. Left panel: reconstructed coupling function A{ip). Right panel: reconstructed first 
derivative of the kinetic function, K'{x) = dK/dx- 


K' must decrease if x increases. Next, Eq.(5.8) gives ~ —62 so that increases with 

|e2|- 

For the model 1 we find that K' slightly decreases from z = 0 to z ~ 4. This will lead 
to a large deviation from the A-CDM matter clustering, as the factor ei defined in Eq.(2.13) 
remains unsuppressed by the denominator K' until 2 ~ 4. This can also be seen in the upper 
curve of the right panel in Fig. 2. However, one cannot take XA too small to achieve low values 
of K' z ^ 2. Indeed, to have a well behaved K-mouflage scenario, without ghosts and small 
scale instabilities, we must have AT' > 0 and K' + 2xK" > 0. Moreover, the reconstruction 
developed in section 5.2 breaks down if K' becomes negative, since a positive function U{a) 
cannot be consistent with a negative K' in Eq.(5.7). These constraints define a lower bound 
on the parameter ja, for given values of the other parameters. 
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6.2.2 Reconstructed coupling and kinetic functions 

The previous reconstruction of the time evolution of cubic models from the {U{a),A{a)} 
parameterisation allows us to rebuild the held dependence of both A{lp) and K'{x)- This is 
shown in Fig. 4. The coupling function A{i~p) is essentially linear over the range of background 
held values (p probed by the models in the past of the Universe. This is consistent with the 
fact that its hrst derivative (5 is almost constant and close to 0.1, as seen in the left panel 
of Fig 2. As explained in section 6.2.1, this is due to the facts that we considered smooth 
functions U{a) and A{a) and that the scalar held excursions are not large, (p/Mpi < 0 . 6 , as 
seen in the right panel in Fig. 1. 

On the other hand, the reconstruction of K'[x) clearly shows that the kinetic function 
is very sensitive to the value of the parameter XA- Because of our common normalization of 
the models at z = 0, x(fl) ^^cid K'{a) converge at low z to common values for all three models, 
as seen in Fig. 3. This implies that K'{x) = 1 for the common value yo — 2 x 10 “^ today. 
At high values of x fo the early matter era, the kinetic function reaches a cubic form, whence 
K' ^ X^. This transition happens at much lower redshifts for the model 3, with a high 7 ^, 
than for the model 1, with a low XA- This is consistent with Fig. 3, where we found that a 
smaller 7 ^, whence a greater |e 2 | at intermediate redshifts, leads to a greater x ~ |c 2 p and 
a lower K' ~ l/|c 2 |- Therefore, as we decrease 7^4 the curve K'{x) shifts to the lower right 
corner in Fig. 4. In particular, the almost constant value K'[a) ~ 1 until z ~ 10 found in the 
right panel in Fig. 3 for model 1 translates into a plateau with AT' ~ 1 until x — 10 in Fig. 4. 
This gives a model where there is no K-mouflage screening of the cosmological background 
until z ~ 10. At very high values of x we can see a change of slope of the kinetic function. 
This is because we enter the radiation era, where the relation between the exponent va of 
62 ( 0 ) and the exponent m of K{x) changes from (5.32) to (5.33). With our choice of a fixed 
VA set by the matter era relation, m = 3 and iza = 3{m — l)/(2m — 1) = 6/5, we hnd in 
the radiation era mrad = 9/8 so that K' ~ x^/® for x > Xeqj where Xeq is the value reached 
at equality between the matter and radiation components. Since the K-mouflage model is 
very close to a A-CDM cosmology at these high redshifts we only reconstructed the kinetic 
function K'(x) up to 2 ; = 10“^ in Fig. 4. 

6.3 Comparison with A-CDM 

The K-mouflage models that we have reconstructed in the previous section show deviations 
from A-CDM at both the background and perturbative levels. We compare the K-mouflage re¬ 
sults to the A-CDM predictions in Figs. 5-6, where all scenarios obey the same normalization 
for the cosmological parameters today. 

We can check in the left panel in Fig. 5 that the cosmological density parameters Dm( 2 ) 
and D(je(-^) remain close to the A-CDM predictions, with percent deviations as for other 
background quantities. This order of magnitude is set by the value of 62 , 0 - In agreement with 
previous works, the dark energy density becomes negative at high redshift in the K-mouflage 
scenarios, which yields Dm > 1 and Dje < 0. However, this only occurs at high z when 
the dark energy density is subdominant. Therefore this feature is not a problem as > 0 
always. This change of sign of pde corresponds to the divergence and change of sign of the 
effective dark energy equation of state parameter w^e = Pde/pde shown in the middle panel 
in Fig. 5. At low z, when the dark energy density becomes important, all models converge 
to Wde — as we remain close to the A-CDM cosmology. Notice that the dark energy 
equation of state can cross the phantom divide w^e < — 1 even though the propagation speed 
cl > 0, as shown in the right panel in Fig. 8 below. This is a drastic difference with iC-essence 
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Figure 5. Left panel: matter and dark energy cosmological density parameters, as a function of 
the redshift. We show the three K-mouflage models of Table 1 and the A-CDM reference cosmology. 
Middle panel: effective dark energy equation of state parameter Wde = Pde/Pde for the K-mouflage 
models. Right panel: relative deviation of the Hubble rate from the A-CDM reference, for the three 
K-mouflage models. 


models [43] due to the coupling to matter of K-mouflage. We can check in the right panel 
in Fig. 5 that the Hubble rates converge to the same normalization Hq today. As we move 
to higher redshift, the K-mouflage Hubble rates first grow more slowly than in the A-CDM 
scenario, which yields a negative deviation Aif at 2 ~ 1. This is due to the decrease of the 
K-mouflage dark energy density with redshift until it vanishes and next becomes negative, 
as explained above, whereas the A-CDM dark energy density remains constant (the matter 
densities remain identical, po/a^, in all scenarios). The K-mouflage Hubble rates become 
greater than the A-CDM expansion rate at high z, when the dark energy density is negligible. 
This is due to the time dependence of the K-mouflage Planck mass (in the Jordan frame). 
Indeed, all Planck masses are normalized to the same value Mpi today, and we have seen 
in Eq.(3.13) that the K-mouflage Planck mass is given by Mpj = /Mpj = A“^Mpp This 
gives a Newton’s constant that grows with redshift and leads to a greater Hubble 

expansion rate from the Friedman equation. The time dependence of the K-mouflage Newton 
constant can be read from the time dependence of the conformal factor A shown in the left 
panel in Fig. 1. For instance, in the model 1 with the greatest deviation from A-CDM the 
Newton constant is higher by 10% at high redshift, z > 100, than its value today. 

In agreement with expectations, the deviations from the A-CDM background quantities 
decrease from model 1 to 2 and from 2 to 3. This corresponds to smaller values of the factors 
ei and 62 and to a higher kinetic function derivative AT', that is, to a stronger K-mouflage 
non-linear screening of the cosmological background. 

From the reconstruction of the K-mouflage model, and more precisely from its Hubble 
expansion rate and the factor ei defined in Eq.(2.13), we can obtain the linear growing and 
decaying modes D±{a) of the matter density held from Eq.(2.14). This also provides the 
linear growth rate /+(a) and the growth rate parameter 7+(a) dehned by 


with /_|_ = 


/+(a) 


din D_|_ 
din a ’ 


7-r(o) 


fo/+ 

InDm ’ 


( 6 . 6 ) 
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Figure 6. Left panel: relative deviation of the linear growing mode D+{z) from the A-CDM 
reference. Middle panel: relative deviation from the A-CDM reference of the linear growth rate 
/_i_ = dlnD^/dlna. Right panel: growth rate parameter ^+{z) for the K-mouflage models and the 
A-CDM reference. 


We show the relative deviations of Dj^[z) and f+{z) from the A-CDM predictions in 
Fig. 6 . The deviation of the linear growing mode Dj^ can reach up to 9% at ^; = 0, even 
though 62 = —0.01 today. This is because the impact of the modified gravity is cumulative 
with time on the linear growing mode. This can yield a significant amplihcation factor, as 
compared with the magnitude of € 2 , if K-mouflage screening is delayed to high redshift, as 
for model 1. The deviations of the growth rate /+ are of the order of a few percents. Again, 
as for the background quantities, the deviations of the perturbed quantities from the A-CDM 
predictions increase from model 3 to 2 and from 2 to 1. This corresponds to higher ei and 62 
and lower K'. In terms of the parameters of Table 1 this corresponds to a lower 7 ^. 

The right panel in Fig. 6 shows that the behaviour of the growth rate parameter 7 +( 2 ) 
obtained in the K-mouflage models is very different from the one found in the A-CDM cos¬ 
mology. This is because of two effects: a) the fifth force amplifies the growth of structures 
and b) the dark energy density actually becomes negative at high redshift in the K-mouflage 
scenarios, which corresponds to Dm > 1 and to a change of sign of In Dm. At low z, as in 
A-CDM scenario, we have /+ < 1 and Dm < 1, and because we always remain close to the 
A-CDM cosmology we have 7 + ~ 0.55. As we go to higher redshift, we converge to the 
Einstein-de Sitter cosmology where both /_|_ and Dm go to unity. Therefore, both ln/+ and 
In Dm become very sensitive to the details of the cosmological model and their sign can change 
with the models. In the A-CDM case we always have /+ < I and Dm < I, whence 7 + > 0, 
and the parameter 7 + does not evolve much. In contrast, in the K-mouflage scenarios we have 
/+ > I at high redshift because of the acceleration of structure formation by the fifth force. 
This leads to a change of sign of 7 +, which becomes negative at 2 ; ~ 2.3 for model I, and at 
higher redshifts for models 2 and 3. Next, at still higher redshift {z > 3 for model I) the dark 
energy density becomes negative while Dm becomes greater than unity. This gives rise to a 
second change of sign of 7 +, which turns positive again, after it diverges by going from —00 
to -|-oo. The behaviour of 7 +(-z) is thus a distinct signature of these K-mouflage scenarios. 
However, there is no dramatic change to the linear growing mode nor to the linear growth 
rate /+ and this only relies on the change of sign of ln/+ and In Dm while /+ and Dm are 
close to unity in the matter dominated era. In practice, it is probably more robust to measure 
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Figure 7. Left panel: effective function A(a) in units of p^. Right panel: effective function c(a), 
multiplied by the factor /po¬ 


und use independently nm(^) and f+(z} in the data analysis, as the exponent 7 + is likely to 
be rather unstable because of this high sensitivity to small deviations in the matter era. We 
leave to future works a detailed analysis of these points. 

6.4 Effective time dependent functions 

The effective action depends on several operators with coupling functions which are time 
dependent. For the case of K-mouflage scenarios they can be expressed in terms of the K- 
mouflage kinetic and coupling functions K(x) and A((f') as in Eqs.(3.7) and (3.13), which 
give 


f{a) = A-\ M^,ia) = A-^Mp,, A{a) = A-^M\xK'- K), 
-f\a)c{a) = A-^M^xK' - = A-^M^x^K". 

i im 

We show these effective functions in Figs 
by [42] 


The function A(a) generalises the cosmological constant to this dynamical setting while 
c(a) represents the effects of the kinetic energy of the K-mouflage field. In Fig. 7, we find 
that the cosmological evolution of A(a) only starts at ~ 10 when dark energy becomes 
dynamical. This corresponds to the time when —g^^(a)c(a) becomes of order unity and K(x) 
starts to deviate from the low-y; constant regime .K ~ — 1 . In the early matter era A(a) and 
—g^^(a)c(a) grow with redshift as K XK' ~ A™- The difference between the behaviours 
of A(a) and —(f^{a)c{a) at low redshifts comes from the fact that A(a) is sensitive to the 
non-zero value of the kinetic function at A = 0 ) which plays the role of the cosmological 
constant, whereas c(a) is only sensitive to the kinetic term K'x and decays as A at late times. 
Again, the deviations from the A-CDM behaviour, a constant A(a), are greater for the model 
1 , which corresponds to smaller XA- 


. 7 and 8 , as well as the propagation speed c^, given 
K' 


K' + 2xK" 


( 6 . 8 ) 
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Figure 8. Left panel: effective function M^ia) multiplied by the factor Right panel: 

scalar propagation speed (in units of the speed of light). Notice that it is always positive. 


At high redshift the quadratic mass term shows the same growth with z as 

A(o) and —g^^(a)c(a), as for a kinetic function that has a power law asymptote at high 
X. For the model 1 this term actually becomes negative at low redshift, with the change of 
sign of K", in agreement with the reconstructed kinetic function displayed in Fig. 4 [we can 
see in that figure that K'{x) has a negative slope over 10“^ ^ A ^ 0-5]- This implies that at 
low redshift the scalar propagation speed Cg is greater than the speed of light for this model, 
as we can check in the right panel in Fig. 8. At high redshift, and at all times for the other 
models where K'{x) is a monotonic increasing function, Cg is smaller than the speed of light. 

This again illustrates the different behaviours obtained for K-essence and K-mouflage 
models. Whereas in K-essence scenarios has to be greater than unity in some redshift 
range in order to reproduce a A-CDM-like expansion history [47], in K-mouflage scenarios 
we usually have < 1 around the cosmological background at all redshifts, except for a 
subclass of models such that K" < 0 at low x which have > 1 at low redshifts. These 
different behaviours come from the fact that the evolution equations of the scalar field are 
different. In the K-essence scenario the Klein-Gordon equation has no external source term 
and shows several fixed points, which also depend on the current expansion era. Then, the 
evolution of the Universe (as it goes from radiation, matter and dark-energy eras) drives the 
scalar field from one fixed point to another [38]. In the K-mouflage scenario, the scalar field 
is driven by the coupling to matter, which gives rise to an external source term in the Klein- 
Gordon equation (2.10). Another important difference is that the dark-energy equation-of- 
state parameter tCde of Eq.(5.15) is different from the scalar-field equation-of-state parameter 
Wip of Eq.(5.16). This leads to different behaviours for Wde{z) and c^{z). 

The superluminality encountered in the case of model 1 in the late-time Universe does not 
entail any causality issue in the present Universe in astrophysical situations where Newton’s 
potential is small enough, see the appendix B for more details. 

7 Conclusion 

We have described how K-mouflage models can be expressed within the framework of the 
effective field theory of dark energy. Because these K-mouflage scenarios are fully defined by 
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their nonlinear Lagrangian, we obtain a fully nonlinear effective-field-theory action that can 
be expanded up to all orders. In the unitary gauge, we find that the K-mouflage scenarios 
only generate the operator at each order n, whereas more general scenarios can also 

generate additional operators involving the extrinsic curvature of the constant time slicing 
surface. This is a distinct signature of K-mouflage models within the class of models coupled 
to matter If observations analysed within the effective held theory framework show that 
such additional operators are required to match data, this would at once rule out all these 
K-mouhage models. 

At the quadratic order, which fully determines the cosmological background and the 
linear regime of the growth of large-scale structures, the effective action is dehned in terms 
of time dependent functions, /(r), A(r), c(r) and M|(r), and we have shown how they can 
be related to the K-mouhage kinetic and coupling functions K{x) and A{(p). We have also 
described how the usual equations of motion associated with K-mouhage models are recovered 
within the effective held theory approach. In particular, at the linear level the scalar held vr 
that appears within this formalism is related to the original K-mouhage scalar held by a 
relation of the form vr oc 5^p. However, the two helds are not identical. Although we obtain 
a closed form for the exact nonlinear equation of motion of the held vr, it has a signihcantly 
more complex form than the one of the original held <p. 

Next, we have presented a simple parameterisation, in terms of two time dependent 
functions U{a) and A{a), that fully determines the cosmological K-mouhage behaviour. One 
advantage of this parameterisation, as compared with the dehnition in terms of K{x) and 
A((p), is that it provides a simple reconstruction of the past cosmological behaviour from 
the normalization of the cosmological parameters today. In contrast, given generic functions 
K(x) S'lid A(y?), because the equations are non-linear one must apply some iterative search 
algorithm to set the mass scale A4^ in order to reach a given dark energy cosmological 
parameter today (which means computing the full background evolution from a <C 1 until 
today for each trial). 

By explicitly choosing the function A(a), and therefore € 2 ( 0 ), we also ensure that some 
cosmological and Solar System constraints are satished. More precisely, le 2 ,o| < 0.01 ensures 
that the current drift of Newton’s constant remains below the upper bound obtained from the 
Lunar Ranging experiment. Moreover, we typically have |e 2 I ~ ei = 2 fi'^/K', see Eqs.(5.35) 
and (2.13). At late times we also have x — ^ ^ind K' ~ 1, which gives |e 2 | ~ ei ~ 2/3^. 
Therefore, |e 2 ,o| < 0.01 also implies /3 < 0.1. This ensures that the upper bound on the 
deviation from the Hubble expansion rate predicted by the A-CDM cosmology at Big Bang 
Nucleosynthesis is satisfied (this deviation is not due to the dark energy component, which is 
subdominant at high redshift, but to the time dependence of the Planck mass). It also means 
that deviations from the A-CDM matter power spectrum are of the order of a few percents 
on large linear scales, which remains consistent with observations. 

By choosing an exponent m > 1 for the highly non-linear regime of the kinetic function, 
K{x) ~ X^, which corresponds to € 2 ( 0 ) ^ 3)/(2m 1 ) matter era within the 

parameterization {U,A}, we also ensure that the dark energy component is subdominant 
at high 2 ; and we recover the Einstein-de Sitter expansion in the matter era. Then, all 
cosmological constraints are satished at these qualitative and semi-quantitative orders. A 
more precise study, through a combined MCMC analysis using various cosmological probes. 


^N-essence would also have the same type of expansion but its features are drastically different, for instance 
the slip parameter 7 (see appendix A.2)) would only differ from one only for K-mouflage. 
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would provide the exact parameter space of K-mouflage models that is allowed in terms of 
the parameters of Table 1. 

However, we must point out that the analysis presented in this paper does not tackle 
all Solar System constraints. In particular, the shape of the kinetic function K(x) at large 
negative x is strongly constrained and must satisfy K' ^ 1 and ^ ■ This is beyond 

the reach of the effective field theory approach, which can only probe the kinetic function 
K(x) in thn cosmological regime y > 0. Thus, the K-mouflage scenarios provide a simple 
example where the cosmological and Solar System regimes are mostly decoupled, although 
there are some common constraints as explained above, such as the value of the coupling 
strengh /3. 

The results presented in this paper can be useful for numerical studies of cosmological 
probes at the background and linear levels. In particular, numerical codes that have been 
designed for the effective field theory framework can use the functions /(r), A(r), c(r) and 
(t) obtained in this paper to compute the observational signatures of K-mouflage models. 
This is left for future work. 
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A Effective-field-theory action at quadratic order 
A.l Scalar dynamics 

The unitary gauge is useful to determine the operators which define the effective action of 
K-mouflage. In particular, the action (3.5) shows that only the operator is generated 
by the K-mouflage models. On the other hand, it hides the presence of the scalar degree of 
freedom which modifies gravity and affects the growth of structure. As recalled in section 3.3, 
the scalar field can be unravelled by the transformation (3.14), which restores the Lorentz 
invariance through the explicit introduction of the scalar field vr. 

In section 3.3 we considered the fully nonlinear longitudinal gauge action (3.11), which 
led to the exact action (3.16). However, in practice one uses effective field theories to describe 
theories defined up to some finite perturbative orders, including all operators that are allowed 
by symmetries up to that order. This has the benefit of providing a very general framework, 
which can encompass many different theories, at the cost of leaving the nonlinear closure 
unspecified. In this spirit, if we only consider the equations of motion at linear order over 
the field fluctuations, we only need the effective action up to quadratic order over the metric 
potentials $ and 'L defined in Eq.(2.12) and over the scalar field vr. Notice that r coincides 
with the time-slicing of the background cosmology. In the small-scale and non-relativistic 
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limits considered in this paper, this gives 


52 =y (fx 
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-34'^ + (V^')^ - 2V^> • 4 ( Svr^' + Vvr • V(^> - 2^ 
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4 

which also reads as 
S2= [ d^xa^lfM^i 


(-TT^ + (V7r)2) 


WfM^X 


-3^ + (V^) - 2V$ • ^ 371-^' + Vtt • V(^> - 2^' 
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+ c [tt^ - (Vtt)^] + - a^$(5p|, 


/ 


{K' + 2xK")n‘^ - K'{V-Kf 
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-a^^dp}, 


f 


(A.l) 


(A.2) 


where 5p is the pressureless matter perturbation in the Jordan frame. In the first expression 
(A.l) we wrote the quadratic action in terms of the original K-mouflage kinetic function K{x), 
whereas in the second expression (A.2) we used the effective-field-theory functions /(r), c(r) 
and which are given by Eqs.(3.13) or (6.7). As expected, in the Jordan-frame action 

(A.l) there is some mixing between the scalar field vr and the gravitational potentials <I> and 
\k. This corresponds to the fact that the metric potentials explicitly depend on the scalar 
field TT in Eq.(4.17). 

To diagonalize the kinetic terms it is convenient to introduce the Einstein-frame poten¬ 
tials, 

f - f 

$ = 4> -h —TT, T = T - —TT, (A.3) 

2/ 2/ ’ ^ ^ 

and the quadratic actions (A.l) and (A.2) become 


52 = / d^x aU /M^i -3T^ + (VT)^ - 2V4> • VT 


+ [{k' + 2|K")^2 - k'iVTrf] - U - ^vr j Jpj, 


(A.4) 


and 


52 = / aU fMk -3T^ + (VT)^ - 2V$ • VT 


/■ 

+ (17^Pi + ^ Jpj, 


(A.5) 


Using this action, we can immediately deduce that the scalar field tt couples to matter with 
a strength f/2f. Moreover the field tt propagates with a speed given by Eq.(6.8) in the 
main text, which also reads in terms of the effective-field-theory functions as 


= 


Mf. A,f2 „ 

4 / ^ 


^M2^ + c + 2a-2M| 


(A.6) 


[Let us recall that the dot denotes a derivative with respect to conformal time r and we 
use comoving spatial coordinates.] This expression coincides with the one obtained in [4] for 
general effective field theories, and with the one obtained in [37] for K-mouflage models (in 
the small-scale and high-frequency limits). 


31 















A.2 Perturbations 


We now consider the dynamics of cosmological perturbations focusing on Cold Dark Matter. 
The matter energy momentum in the Jordan frame is conserved, 


= 0 , 


(A.7) 


implying that the continuity and the Euler equations coincide with the ones in GR. The 
modifications to the dynamics of perturbations are solely induced by the change to the Pois¬ 
son equation for the Jordan-frame metric potential We focus on the sub-horizon regime 
of perturbations where in terms of Fourier modes k/a S> H, implying that the time varia¬ 
tions of TT and can be neglected compared to their spatial variations. In this quasi-static 
approximation, the effective action (A.5) reduces to 


(V7r)2 - ^l> - 5p 

(A.8) 

where we have expressed all factors in terms of the effective-field-theory functions /(r) and 
c(r). The resulting Poisson equation for the Einstein-frame Newtonian potential is not mod¬ 
ified and we recover from the action (A.8) 


S2= I (fx a^\ fMly (VT)2 - 2V^ • VT 


4/ 


AJpj -|- c 


- 

$ = T = Tn with —tj-'I'n 


2M2j’ 


(A.9) 


where Mpj = fMp^ is the Jordan-frame Planck mass and 6p the Jordan-frame matter density 
fluctuation. The field vr satisfies a similar relation. 




-TT = —- 


4/ 


iff ^ -1- c 
4/ ^ 


bp 


(A.IO) 


As a result the Poisson equation for the Jordan-frame metric potential $ is modified and 
becomes 

^Mpi + c Sp 




-$ = 


+ c 


(ATI) 


4/ "'"PI 

which also writes in term of Newton’s potential, Tn = 'I'N) as 


jrAJp, 2/32 

$ = (l-bei)TN with ei = -= 

and we recover Eq.(2.13). The effect of this modification of the Poisson equation appears 
naturally in the growth of the matter density contrast b = bp/p, which satisfies in the linear 
regime 

0 + 2H^ - + ei)b = 0, (A.13) 

in agreement with Eq.(2.14). This means that the growth of large-scale structures defines a 
modified Newton constant as 


11 TVf 2 _|_ g 

= 0/2 ~ 0 -0N = (1 + ei) ^N, (AT4) 

+ c 
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where = I/SttM^j = l/87r/(r)Mp[. 

Lensing effects are not sensitive to ‘h but to + 'I' = 2'1'N) implying that the leasing 
potential feels the modified Newton constant 




(A.15) 


ft is one of the main features of these models that gravity does not have a unique signature 
but manifests itself differently for matter and photons. 

Similarly, the ratio between the two Jordan-frame Newtonian potentials 'I' and defined 
by 7 = ^, is not equal to one but 


d' _ 1 - ei 
~ 1 + ei 


1 _ 

K' 

l + Wl' 
r -r 


(A.16) 


2 m 

Again we find a deviation from unity which depends on These results coincide with the 
ones deduced in the main text. 


B Superluminality and Causality 

In this appendix, we comment on the superluminality of scalar perturbations and its link with 
causality. This issue was analysed for K-essence models in [44]. We follow a similar method 
here. Let us first expand the K-mouflage action to second order in cr, where ip = (p + a. Here 
(p{x, t) is a background configuration that may depend on scale and time. We also denote 
X = —g^''d^(pdy(p/2Ai'^^ K' = K'(x) and K" = K"{x)- The second order part of the action 


reads 


52 = ^ / (fxy/^ -K'd^^ad^a + ^{d^^ipd^af 

(B.l) 

ft is convenient to define the disformal metric 

K” \ 


Gf^^ = 7-1 ( K’g^^'^ - 


(B.2) 

with 



7 = + 2xK'' f/‘^ > 0. 

(B.3) 


We consider models for which K' > 0 and K' -|- 2xK" > 0 over all y, so that the metric 
is well defined, see also (5.4). Defining the inverse matrix and the determinant 
G = det(G^j/), we have 

G = 7"5, (B.4) 

where g = det[g^y), and the second-order action can also be written as 

52 = - i/ d^xV^G>^^d^ad^a. (B.5) 

The disformal metric G^u is the metric felt by the scalar perturbation and it is Lorentzian. 
Initial-value problems for a are well posed on any smooth space-like Cauchy surface S for 
the metric G^u, and the solution is unique and propagates causally (see Sec.10 in [48]). 
Around the cosmological background we recover the hyperbolic Klein-Gordon equation with 
the propagation speed Cg of Eq.( 6 . 8 ). 
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Space-times equipped with the metric are stably causal, i.e. there are no time-like 
closed loops, provided there exists a globally defined function / on all space-time which is 
time-like, i.e. d^fd^f < 0 [48]. Following [44], we search for a “global time” / that applies 
to both geometries and G^^, and thus implies the absence of closed causal loops. A simple 
choice is the cosmic time t, which clearly satisfies the required property for the metric 
Considering the Newtonian gauge, which describes all systems where gravity is weak enough 

ds^ = g/j.ydx^dx'^ = —(1 -|- — 2'^T<i)dx‘^, (B.6) 


where 'I'n is the Newtonian potential, we have 

g^'^d^tOyt = ^ < 0 for 4 'n > -1/2. (B.7) 

i -h 

Since we focus on systems with |'kN| ^ 1 (e.g., ~ 10“® in the Solar System), we have 

g^’^d^tdyt < 0. On the other hand, we obtain 


G^^’^d^tdyt = 


K'jl + 2 ^ n ) + K''ido^)yM^ 

7(1 + 24'n)2 


(B.8) 


and therefore ^ 

G^^^d^tdyt <0 for C = K' + > 0, (B.9) 

where we used the approximation 1 -|- 2'1 'n — 1- Around the cosmological background, where 
X = {d(p/dt)"^/2M^ , we obtain C = K' + 2xK", hence C > 0. Around a static background, 
we obtain C = K' whence C > 0. 

For more general backgrounds, we can see from Eq.(B.9) that C > 0 as soon as K" > 0. 
When K” < 0, we can only conclude in some specific cases. First of all consider cosmological 
situations where 3> (V(^)^, we have C ^ K' + 2xK" > 0. Similarly in almost static 
configurations where {(p^ ^ we have C > K' + K" ^ K' — 2xK". In this 

almost static limit, where X < 0 and is a large negative number, we must apply the restriction 
obtained from the bound on the perihelion advance of the moon, i.e. we require that | Al^^xl ^ 
K' when X ^ 1 [35]. This implies that C > 0 for the models passing the solar system tests. 
We now examine the intermediate cases when x r\j — {d(pjdt)'^l2M-^ ~ —(V(/?)^/2A4'^, i.e. the 
system is both spatially and time varying. Then, we have C ^ K' + 2K"x > 0. In conclusion, 
even when K” < 0, these three cases cover all of the astrophysical situations of interest, and 
a violation of causality is not present. 

Thus, we conclude that g^'^d^tdyt < 0 and G^'^d^idyt < 0 and there are no closed causal 
loops around usual astrophysical and cosmological backgrounds with 'I'n > —1/2. This anal¬ 
ysis fails close to neutron stars or black holes, where 'I'n becomes large, but this is not related 
to the K-mouflage model as it already appears in the metric g^y. More generally, in such theo¬ 
ries where super luminal velocities can be found, it may be possible to find inhomogeneous and 
time-dependent backgrounds, such as boosted bubbles moving apart [46], where closed causal 
loops exist. Moreover, global properties of space-time may break causality. However, these 
two issues also appear within General Relativity itself and one has to resort to the chronol¬ 
ogy protection conjecture, which assumes that the laws of physics prevent the formation of 
backgrounds leading to closed causal loops. Then, as argued in [44, 49], from the point of 
view of causality theories such as K-essence, or K-mouflage in our case, are similar to General 
Relativity and do not lead to greater causal paradoxes. From the perspective of this paper, 
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where we restrict ourselves to linear perturbations around the cosmological background, we 
have shown from the analysis above that causality is preserved, even in the case of model 1 
where > 1 at low redshifts. Within the effective field theory approach, we only reconstruct 
the regime X > 0; the regime y < 0, where problematic backgrounds may appear, is not 
reconstructed and beyond the scope of such frameworks. 
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